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ABSTRACT 

Wilson’s GF matrix method is generally utilized to 
obtain normal vibrations and normal coordinates of 
molecules. The function of the product of G and F in 
Wilson’s GF matrix method is also a key to determine 
the line-shape function (LSF) which expresses 
molecular absorption or emission spectra. In this 
paper, the method for calculating the function of the 
product of G and F is shown. 

Keywords: Line-shape function, GF matrix method, a 
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I. INTRODUCTION 

The line shape-function (LSF) plays an important role 
in a variety of research (ordinary optical absorption or 
emission, band shape of the circular dichroism, 
excitation energy transfer between ions in solids, 
electron transfer between molecules in the solvent and 
in photo synthetic systems, intersystem crossing and 
internal conversion in a molecule, thermal ionization 
of an electron trapped on an impurity in the solid, 
intersystem crossing process of organic light-emitting 
diodes and so on) [1-17], since these transition rates 
by the familiar Fermi golden rule are expressed in 
terms of LSF. If an initial state of system exists in a 
thermalized condition, the LSF is defined as 
following: 

F{co,P)^JJ^p iv {P)Sf ev S 

U V 




-CO 


( 1 ) 


distribution of initial vibrational states, a vibrational 
energy level measured from the bottom of the 
adiabatic potential surface in each electronic state, the 
energy gap between the two bottoms, and the square 
of the vibrational overlap integral between one 
electronic state and another electronic state called the 
Franck-Condon factor, respectively. And the 
independent variables of the LSF, a> and /?, mean the 
angular frequency of the light emitted (absorbed) 
between two energy levels and the value \l(k^T), 
respectively. Here k§ is the Boltzmann constant and T 
is temperature. Namely, the LSF is a function 
depending on temperature and the energy gap 
between two vibronic states. 

Another expression of LSF was derived on the basis 
of the Liouville-von Neumann equation for the 
density matrix [18], and an approximate expression 
was also described. One feature of the approximate 
formula is to express the LSF in terms of both 
changes in molecular structure and in force constants 
between the initial and final electronic states. The 
change in molecular structure is expressed as the 
change in G matrix, and the changes in force 
constants are expressed as the change in F matrix 
between the initial and final electronic states. Here G 
and F are defined in Wilson’s GF matrix method [19]. 
Another feature is that the outline of spectrum shape 
can easily be estimated if both G and F are obtained 
from quantum chemical calculation. 


Here the suffix i or j and the suffix u or v are used for 
specification of a zero-order electronic state and a 
vibrational state, respectively. The variable p w (/f, E iu 
(or Ej V ), A E e , and S iuj f represent the Boltzma nn 


However, in the approximate formula there is the 
following function of the product of G and F: 
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S i =B i AX i . 


(3) 


where Bj is the (3N-6)*3N matrix whose elements 
depend on the molecular structure. Here, i represents 
the i-th electronic state, and N is the number of atoms. 
The G matrix is defined by G, = Zf Af 1 B, J (M is the 
3N*3N diagonal matrix whose elements are the 
nuclear masses). The F t is a symmetric matrix 
consisting of force constants. 

To describe the normal vibration, we employ the 
transformation from the internal displacement to the 
normal coordinates: 


Si — Li Qi 


(4) 


The transformation matrix L, is then determined so 
that it can satisfy the following eigenvalue equation 
and orthonormal condition: 


G, Ft L t - L i A iy 

L? Gf l Lt = I. 


(5) 

( 6 ) 


Here A t is the diagonal matrix whose diagonal 
elements consist of the square of the angular 
frequencies {cy } of the molecular vibration. The 
method for obtaining the matrix L, concretely is given 
in [19-21], 

III. FORMULA OF LSF WITHIN A GAUSSIAN 
APPROXIMATION 

In [18], the approximate formula of LSF within 
Gaussian approximation was derived using Wilson’s 
GF matrix as follows: 


h 2 


Here the suffix e means the electronic excited state. In 
this paper, the calculating method for (2) is 
introduced. 

II. THE DEFINITION OF WILSON’S G- 
MATRIX AND F-MATRIX 

In this section, the definition of the G and F matrix of 
Wilson’s GF matrix method [19] is shown. Within a 
small displacement, the internal displacement 
coordinates S can be related by linear transformation 
of the Cartesian displacement coordinates: 


F(a>,fi)= -—eap 
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(16) 


Here, Aff eq is the structural difference between the 
initial and final state in the internal coordinate. 

IV. A METHOD TO CALCULATE FUNCTIONS 
OF THE PRODUCT OF G AND F 

This section is the main part of this paper and the 
calculation method for functions of the product of G 
matrix and F matrix. Concretely, the solving method 
for (16) is elucidated. 


From Eq. (5), we obtain 
GF= L A U l . 


(17) 


Here, the suffix representing the electronic state is 
omitted. Consequently, we also obtain 
(GFf = LA*L~\ (18) 

Next, we consider the polynomial expansion of 
AGF ): 

f[GF) = CoI+C\{GF)+C2(GF) 2 +C 3 (GFf+- (19) 


From (18), (19) can be expressed as 

f(GF) = L(f j C /J A M )lT l = Lf(A)L~ l . 

M =0 


( 20 ) 
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Because the matrix A is a diagonal matrix, 

[/(/1)],j-/(®, 2 )4, (21) 

Here Sy is Kronecker delta. From Mclaurin series 
expansion for exponential functions, the function 
“coth” can be calculated using (21). 

Then, because the derivative of Jx does not have a 
defined value at 0, Mclaurin series expansion for -Jx 
may immediately seems to be impossible. The case of 
Vgf is the same as the case of 4x . 


However, we consider the polynomial expansion for 
yfx +1 instead of 4gf . Here / is the unit matrix. 
Namely, 

4x +/= c 0 i+c t (x)+c 2 ( xf + (22) 

Then, squaring both side of (22), we obtain the 
following: 

X +1 = C 2 / + 2 C a C x X + (C 2 + 2 C 0 C 2 )X 2 ^23) 

+ 2(C 0 C 3 + C l C 2 )X 3 + ■■■. 

From this equation, it is found that Co=l, Cr=l/2, 

C 2 = - (1/8), 03=1/16, and so on. Thus, 

'Tx+i (24) 

= / + (l/2)X-(l/8)X 2 +(1/16)0 + •••. 


Here, suppose X to be GF - /, we can make the 
following expression from (24) 

VGF = I + (1 / 2)(GF -1) - (1 / 8 )(GF - 1) 2 (25) 

+{\nb){GF-if + ■■■. 


The right hand side of this equation is just the infinite 
power series for the matrix GF. Namely, 4gf can be 
expressed as the power series and be calculated using 
( 21 ). 


V. CONCRETE CALCULATION 

We calculated the LSF of the one-photon emission of 
SO 2 within Gaussian approximation, using (21). As 
for the temperature T, 300K was used. 



EnusOTi Wave Length fnnrO 

Fig. 1 The LSF of SO 2 based on the formula (7). The 
horizontal axis represents the emission wave length of 
S0 2 , and the vertical axis represents the value of the 
LSF. The peak value of LSF is about 5x10 ” 16 seconds 
at 240nm. 

First, we calculated the electronic ground state and the 
first excited state using ab-initio quantum chemical 
calculation software Gaussian 09 [21], The 6- 
31G(d,p) basis set was used in the calculation of both 
the ground and the excited states. We optimized the 
structure of the ground state with HF and of the 
excited state with CIS (Nstates=5), respectively. 

Second, the G matrix, the F matrix, and the L matrix 
in (5) in both the ground and the excited states were 
obtained according to the calculation method in [22], 

Finally, the LSF was calculated using the calculation 
method described in IV. The result of the calculation 
was shown in Fig. 1. The calculated LSF additionally 
agrees with the experimental data [23], 

V. CONCLUSION 

The solving method of (16) are elucidated, and it was 
found that not only (16) but also functions of the 
product of G and F shown in Wilson’s GF matrix 
method can generally be calculated using the angular 
frequencies which are the eigen values of the product 
of G and F. 
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